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An earlier measurement on the 4 + to 2 + radiative transition in 8 Be provided the first electro- 
magnetic signature of its dumbbell-like shape. However, the large uncertainty in the measured cross 
section does not allow a stringent test of nuclear structure models. The present paper reports a more 
elaborate and precise measurement for this transition, via the radiative capture in the 4 He+ 4 He re- 
action, improving the accuracy by about a factor of three. The ab initio calculations of the radiative 
transition strength with improved three-nucleon forces are also presented. The experimental results 
are compared with the predictions of the alpha cluster model and ab initio calculations. 

PACS numbers: 21.60.De, 23.20.Js, 24.30.Gd, 25.55.-e, 27.20.+n 



The nucleus 8 Be is a classic example of the occurrence 
of alpha clustering [l| in nuclei. Its formation from two 
alpha particles provides an intermediate step in the syn- 
thesis of 12 C [2] from the fusion of three alpha particles 
inside the stars. The nucleus is also the stepping stone to 
understand alpha-clustering in heavier self-conjugate 4n 
nuclei. The dumbbell-shaped nucleus exhibits rotational 
states manifested as resonances in the alpha-alpha scat- 
tering system. The electromagnetic transition between 
the excited resonant states in 8 Be, with spin-parities of 
4+ and 2 + , was reported earlier Q in order to provide 
a test for its alpha cluster structure. The measurements 
were made at two beam energies, on and off the 4 + res- 
onance, by detecting the transition gamma rays in co- 
incidence with the two alpha particles arising from the 
decay of the 2 + final state. However, the measured cross 
section (with an uncertainty of ^33%) and the inferred 
reduced electromagnetic transition rate were not precise 
enough to provide a stringent test for various models like 
the cluster model [J] and ab initio quantum Monte Carlo 
model [5J . The uncertainty arose mainly due to the large 
background of 4.44 MeV gamma rays originating from 
the interaction of the incident beam with the window of 
the chamber holding the helium gas target. The present 
work, using essentially the same method, is aimed at a 
more accurate measurement and also at more beam en- 
ergies straddling the 4 + resonance. The essential aspects 
in this improved measurement are a better pixelisation 
of the alpha particle detectors, a more efficient and seg- 
mented gamma ray detector and a better shielding of 
the gamma rays from the beam-window interaction men- 
tioned above. 

The experiment was carried out using beams of 4 He 
from the BARC-TIFR Pelletron Linac Facility at TIFR, 




FIG. 1: (Color online) Schematic of experimental setup. 



Mumbai at energies of 19—29 MeV. The beam current 
was about 1 pnA on the target. The schematic of the 
experimental setup is shown in Fig. 1. The 7-rays were 
detected in a BGO detector array with a photopeak effi- 
ciency of about 23% at E 7 =8 MeV. The array consisted of 
38 hexagonal cross section detectors, of length 76 mm and 
a face to face distance of 56 and 58 mm (in two groups), 
encased in thin aluminum housing. These were mounted 
in close packed groups of 19 each placed at ~70 mm above 
and below the target. Alpha particles were detected in a 
500 fim thick, annular, and double sided silicon strip de- 
tector (SiSD), with 2x16 8 rings (in left and right halves) 
and 16 sectors @ with separate readouts. The active 
portion of the detector had an inner diameter of 48 mm 
and an outer diameter of 96 mm. The gap between the 
adjacent rings was 0.1 mm while that between adjacent 
sectors was 0.2 mm. The left and right halves of the 9 



2 





4000 


1 1 

0!( 12 C) 


1 




3000 














Counl 


2000 
1000 


- 

<X!( 16 0) 

• 


- 

a ( 12 C) 

•;ao( 16 6) 












14 16 18 20 

E a (MeV) 

FIG. 2: Alpha particle spectra from the scattering on the 
mylar target at Ebeam=22.4 MeV for a particular ring of the 
SiSD. Elastic and inelastic peaks (from the first excited states) 
from the 12 C and le O targets are indicated. 



side were separated by 0.4 mm. 

A chamber was designed to mount the strip detector at 
^70 mm from its centre in the forward direction and to 
hold the target helium gas (purity >99.9%) at ^0.8 bar 
pressure. The gas was isolated from the beam line vac- 
uum using ~1 mg/cm 2 Kapton foils at the entry and 
the exit. The helium gas was filled using a gas handling 
setup similar to that used earlier |3|. Conical heavymet 
shields surrounded the Kapton windows in order to shield 
the BGO detector array from the copious 4.44 MeV 7- 
rays produced in the excitation of 12 C in the windows. 
The chamber had the provision to mount a ladder for 
holding a thick aluminum aperture plate with a hole of 
24 mm diameter as well as thin mylar and carbon foils. 
The aperture plate, when placed at the centre of the tar- 
get chamber, shielded the a-particles scattered from the 
Kapton entrance window and limited the effective beam- 
target interaction zone seen by the SiSD. The aperture 
diameter and the SiSD distance were decided on the basis 
of a Monte Carlo simulation 0] to get a reasonable effi- 
ciency for the detection of two a-particles following the 
radiative capture and subsequent decay of the final state 
in 8 Be. The typical effective target length was about 
20 mm and the efficiency for the 2-a detection from the 
final state was about 35%, after including the effect of 
the various dead zones in the SiSD. 

The energy and timing signals of the SiSD were gener- 
ated from each of the 32 6>-rings (divided into two groups 
of the left and right halves) and 16 ^-sectors. The energy 
signals were sent to voltage sensitive analog-to-digital 
converters. Each timing signal was fanned out into two 
paths one being used to generate the overall particle 
event trigger for left rings, right rings and the sectors 
using a logical OR condition among the corresponding 
signals. In the other path the signals were fed to time to 
digital converters (TDCs) for measuring timing with re- 
spect to the 7-ray detector array. The anode signal from 
the photomultiplier of each of the 38 BGO detectors was 



also fanned out for measuring the energy deposited by a 
charge-to-digital converter and for the timing measure- 
ment with respect to the SiSD using TDCs. A logical 
OR condition among the 38 signals produced the 7-ray 
event trigger. The grand event trigger was generated by 
requiring a fast coincidence between event triggers from 
the left and right halves of the SiSD 0-rings and that 
from the BGO-detector array. The data were collected 
in an event by event mode using a CAMAC based data 
acquisition (DAQ) system Q. A 10 Hz pulser signal was 
fanned out and given to the test input of the three SiSD 
preamplifiers for estimating the dead time of the DAQ 
system. 

The energy calibration of the SiSD detector was done 
using elastic and inelastic scattering of a-particles on 12 C 
and 16 using thin carbon and mylar targets. A typi- 
cal a-particle energy spectrum is shown in Fig. 2. The 
energy calibration was performed over the 256 (16x16) 
6 - 4> pixels. The 4.44 MeV and 6.13 MeV 7-rays from 
excited states in 12 C and 16 populated through the in- 
elastic a-particle scattering were used to calibrate the 
BGO detectors. These measurements were made period- 
ically throughout the experiment in order to track the 
possible change in the calibrations of the a-particle and 
7-ray detectors. A stability within ~1% was witnessed 
over the period of the experiment. 

The data were collected at four beam energies of 19.2, 
22.4, 24.7, 28.9 MeV, spanning the 4+ resonance in 8 Be, 
for the integrated beam charges of 81, 90, 125 and 58 
pnC, respectively. The data was analyzed to extract the 
events corresponding to the 7-ray transition to the 2 + - 
final state in 8 Be and the subsequent 2-a decay of the 
final state. The first condition imposed was the prompt 
concidence among the 7-ray detector, at least one of the 
left rings and at least one of the right rings. This en- 
sured a prompt coincidence also between the left and 
the right halves of the SiSD. The sector timing was also 
demanded to be in prompt coincidence with the 7-ray 
detector with two opposite sectors being simultaneously 
in coincidence. These conditions emphasized on the re- 
quired events because the two a-particles from the decay 
of the final state are emitted at the azimuthal angles dif- 
fering by ~ 180°, neglecting the small momentum kick 
due to the transition 7-ray. Fig. 3 shows an example 
of the time spectrum between the BGO detector array 
and the SiSD at a beam energy of 22.4 MeV showing the 
prompt time peak. The hit-multiplicities of the left and 
the right rings were constrained to one for each and the 
energy deposited in the left and right halves (E^ and E#) 
were constructed from the energy calibrations of the cor- 
responding rings. For getting the 7-ray energy E 7 , the 
BGO detector with the highest 7-ray energy was taken to 
be the primary detector. The energies deposited in the 
neighbouring detectors, which were also in prompt co- 
incidence and contained the leaked shower energy, were 
added to that of the primary detector for each event. An 
event by event reconstruction of the total a-particle en- 
ergy, E sum = E^ + Fir was made with conditions of (a) 
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FIG. 3: Example of the time spectrum between the BGO 
detector array and the SiSD detector rings at Eb oa m=22.4 
MeV. 
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both El and being within a lower and a upper limit 
(~1— 13 MeV) and (b) the reconstructed total energy and 
the total momentum of the two a-particles being within 
a proper two dimensional gate. These conditions were 
guided by the Monte Carlo simulations described below. 

Fig. 4 shows a two dimensional plot of E SMm vs E 7 at 
the beam energy of 22.4 MeV. A band of events around 
an E sum of 13 MeV and E 7 of 8 MeV can be clearly iden- 
tified. These events arise from the radiative capture of 
the two a-particles to the 2 + resonance in 8 Be. A one di- 
mensional spectrum of E tot = E sum + E 7 is generated by 
putting one dimensional gates on E sum of 8.8—15.0 MeV 
and E 7 of 3.4—10.5 MeV as suggested by simulation re- 
sults. Similar E tot spectra were generated at other beam 
energies by putting appropriate gates on these quanti- 
ties. Fig. 5 shows the E tot spectra at all the four beam 
energies. The peaks in the spectra (not apparent at the 
highest beam energy which is beyond the extent of the 
4 + -resonance) , corresponding to the 7-ray transition be- 
tween the resonances, were used in the calculation of the 
capture cross sections. 

The extraction of radiative capture cross sections re- 
quires a simulation of the experimental set up using a 
Monte Carlo code. Such a simulation was done in two 
parts - one for the detection of the two a-particles and 
the other for the response of the BGO array to the in- 
cident 7-rays. The simulation for a-particle detection 
took into account the extended gas target, the aperture, 
the angular distribution of the a-particles emitted after 
the 7-decay and the geometry of the SiSD. The energy 
losses of the beam and decay a-particles, and the an- 
gular straggling, were calculated for each event using the 
SRIM code . The inefficiencies of the SiSD due to gaps 
between rings and sectors were included. The 7-ray re- 
sponse was simulated using GEANT3 [l(| with the angu- 
lar distribution from the aligned 4 + to the 2 + final state 
included. For each event, the 7-ray energy was Doppler 
corrected. The simulated event by event data was writ- 
ten in a file for analysis by the same program that was 
used to sort the actual data. The simulated data were 



FIG. 4: Two dimensional spectrum between E sum =El+Er 
(see text) and E 7 , generated with proper cuts as discussed in 
the text, at E bcam =22.4 MeV. 

sorted to create the E tot spectrum with the same condi- 
tions as used in the case of the actual data. Starting from 
the No events corresponding to the 4 + to 2 + transition 
and the subsequent 2-a decay, the counts N in the same 
peak regions as shown in Fig. 5 were calculated from the 
simulated spectra to get the overall detection efficiency 
(N/Nq) of the experimental set up. The simulation also 
provided the effective target thickness. The capture cross 
sections were extracted using the integrated beam charge, 
the target thickness and the detection efficiency. The ef- 
fective a-particle energy <E Q > at each beam energy was 
also extracted from the simulation after knowing the in- 
teraction region and the energy loss of the incident beam 
in the entrance window and in the target gas up to the 
interaction region. The spread in the effective energy 
due to the finite extent of the interaction region was less 
than 0.14 MeV. The extracted cross sections at the four 
beam energies and the corresponding effective a-particle 
energies are tabulated in Table I. 



TABLE I: Effective a-particle energies <E a > and the radia- 
tive capture cross sections <r 7 extracted from the data at four 
beam energies (E beam ). 
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The extracted cross sections are plotted against the ef- 
fective a-particle energies in Fig. 6. The cross section at 
the resonance energy is consistent with the earlier mea- 
surement Q, but with an error of ~10% as compared 
to the earlier 33%. Fig. 6 also shows the calculated cross 
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FIG. 5: One dimensional spectra of Etot =E sllm + E 7 gener 
ated with proper cuts at four beam energies. 



sections from the cluster model calculation of 4]. The 
contribution from the partial waves of 1=0, 2, 4 are added 
incoherently in this plot. The comparison with experi- 
mental data is good in the rising part of the cross section 
profile but deviates at higher energies. Whether a differ- 
ent choice of the a — a potential along with a coherent 
summing over the various partial waves will improve the 
comparison remains to be seen. It may be mentioned 
that there is some ambiguity in the choice of the poten- 
tials giving similar values for the energies and widths of 
the resonant states of 8 Be. 

Ah initio calculations of the radiative transition 
strengths in 8 Be, using realistic two- and three-nucleon 
interactions, were first reported in These varia- 

tional Monte Carlo (VMC) calculations of the electric 
quadrupole moment Q and B{E2) values indicated that 
the low-lying spectrum of 8 Be is well-described by the ro- 
tation of a common deformed two-a structure. Recently 
it has become possible to evaluate electroweak transitions 
with the more accurate Green's function Monte Carlo 
(GFMC) method while improvements in the three- 
nucleon forces have also been made. We report here new 
GFMC calculations using the Argonne wis (AV18) two- 
nucleon [l2j and Illinois-7 (IL7) three-nucleon |l3| poten- 
tials, which give a very nice reproduction of the energy 
spectrum and other properties of light nuclei in this mass 



200 



150 



tTlOO 



50 



Present data 
O Old data: Ref. 3 
Theory: Ref. 4 




10 15 20 

E a (MeV) 



25 



30 



FIG. 6: (Color online) Extracted capture cross sections plot- 
ted against the effective a-particle energy (see text). The last 
point indicates the upper bound of the cross section. The 
continuous line shows the result of a model calculation. 
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FIG. 7: (Color online) GFMC propagation with imaginary 
time t of the (a) energy, (b) point proton radius, and (c) E2 
matrix element; open stars denote the values extracted from 
the calculation. 



range flij ]. 

An initial VMC calculation is made to generate a start- 
ing wave function, which the GFMC method then sys- 
tematically improves upon by a propagation in imaginary 
time r. The 8 Be 2 + and 4 + excited states are particularly 
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challenging because they tend to break up into two sepa- 
rate a-particles as r increases. Figure El shows the propa- 
gation with imaginary time of the energies, radii, and E2 
matrix elements. In Fig. Eta), the energies of the states 
are seen to drop rapidly from the initial VMC energies 
at r = 0. The + ground state energy stabilizes and is 
well-fit by a constant averaged over r = 0.1 — 0.3 MeV -1 . 
The 2 + state shows a very subtle decrease over the same 
range, while the 4 + state drifts significantly lower; the 
energies quoted below are obtained from a linear fit us- 
ing the value at r = 0.1 MeV -1 , with the Monte Carlo 
statistical error augmented by the range of values from 
t = 0.08 - 0.12 MeV -1 . This choice of r should encom- 
pass the bulk of the improvement in the wave functions 
provided by the GFMC algorithm, before the tendency 
to dissolve into two a-particles sets in. 

This tendency to dissolution is seen more strongly in 
the evolution of the point proton radii shown in Fig.Etb). 
The + ground state radius is flat as a function of r, while 
the 2 + and 4 + states both increase steadily from about 
t = 0.1 MeV -1 . The associated electric quadrupole mo- 
ments, which are not shown, also increase steadily. Fi- 
nally, the E2 matrix elements, shown in Fig. Etc), also 
increase with r, the effect being particularly pronounced 
with the (4 + — > 2 + ) transition. 

Results for the energies E, point proton radii r pi 
electric quadrupole moments Q, and B{E2) transition 
strengths are given in Table [TTJ The absolute energies of 
the states are in excellent agreement with experiment. 
The quadrupole moments and B{E2) values are consis- 
tent with an intrinsic quadrupole moment Qq of 32 ± 1 
fm 2 , which is about ~ 20% bigger than the original VMC 
calculation of @. 



TABLE II: GFMC results 



r 


E [MeV] 


r p [fm] 


Q [fm 2 ] 


B(E2 I) 


0+ 


-56.3(2) 


2.40 







2+ 


-53.1(1) 


2.45(1) 


-9.1(2) 


20.0(8) 


4+ 


-45.1(2) 


2.48(2) 


-12.0(3) 


27.2(15) 



A comparison with the ab initio calculation needs the 
present experimental result to be expressed in terms of 
the B(E2) value for the 4 + to 2 + transition. Whereas 
this is not straightforward, an approximate value can be 
calculated assuming a Breit Wigner form factor for the 
4 + resonance and using the experimental cross section at 
the resonance energy. This gives a partial gamma width 
r 7 =(0.48±0.05) eV and a B(E2) value of 21±2.3 e 2 /m 4 . 
This is somewhat lower than the calculated value. How- 
ever, a better comparison will be possible after perform- 
ing the ab initio calculation as a function of the alpha 
particle energy. The present experimental results, be- 
sides putting the a-cluster structure of 8 Be on a firmer 
footing, will provide data for testing the future calcula- 
tions incorporating the reaction and the structure aspects 
in a seamless manner. 
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